Cytochrome P450-specific human PBPK/PD models for the organophosphorus pesticides: chlorpyrifos and parathion.

Organophosphorus pesticides (OPs) remain a potential concern to human health because of their continuing use worldwide. Phosphororthioate OPs like chlorpyrifos and parathion are directly activated and detoxified by various cytochrome P450s (CYPs), with the primary CYPs involved being CYP2B6 and CYP2C19. The goal of the current study was to convert a previously reported human pharmacokinetic and pharmacodynamic (PBPK/PD) model for chlorpyrifos, that used chlorpyrifos metabolism parameters from rat liver, into a human CYP based/age-specific model using recombinant human CYP kinetic parameters (V(max), K(m)), hepatic CYP content and plasma binding measurements to estimate new values for acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) inhibition and to use the model as a template for the development of a comparable parathion PBPK/PD model. The human CYP based/age-specific PBPK/PD models were used to simulate single oral exposures of adults (19 year old) and infants (1 year) to chlorpyrifos (10,000, 1000 and 100 µg/kg) or parathion (100, 25 and 5 µg/kg). Model simulations showed that there is an age dependency in the amount of blood cholinesterase inhibition observed, however additional age-dependent data are needed to further optimize age-specific human PBPK/PD modeling for these OP compounds. PBPK/PD model simulations estimated that a 4-fold increase or decrease in relative CYP2B6 and CYP2C19 content would produce a 9-22% inhibition in blood AChE activity following exposure of an adult to chlorpyrifos (1000 µg/kg). Similar model simulation produced an 18-22% inhibition in blood AChE activity following exposure of an adult to parathion (25 µg/kg). Individuals with greater CYP2B6 content and lower CYP2C19 content were predicted to be most sensitive to both OPs. Changes in hepatic CYP2B6 and CYP2C19 content had more of an influence on cholinesterase inhibition for exposures to chlorpyrifos than parathion, which agrees with previously reported literature that these CYPs are more reaction biased for desulfurization (activation) and dearylation (detoxification) of chlorpyrifos compared to parathion. The data presented here illustrate how PBPK/PD models with human enzyme-specific parameters can assist ongoing risk assessment efforts and aid in the identification of sensitive individuals and populations.

Human hepatic cytochrome p450-specific metabolism of parathion and chlorpyrifos.

Organophosphorus pesticides (OPs) remain a potential concern to human health because of their continuing worldwide use. Thiophosphorus OPs, once bioactivated by cytochromes P450 (P450s), form oxon metabolites, which are potent acetylcholinesterase inhibitors. This study investigated the rate of desulfation (activation) and dearylation (detoxification) of parathion and chlorpyrifos in human liver microsomes. In addition, recombinant human P450s were used to quantify, for the first time, the P450-specific kinetic variables (K(m) and V(max)) for each compound for future use in refining human physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) models of OP exposure. CYP1A2, 2B6, 2C9, 2C19, 3A4, 3A5, and 3A7 were found to be active to a widely varying degree in parathion metabolism, whereas all, with the exception of CYP2C9, were also found to be active in chlorpyrifos metabolism. CYP2B6 and CYP2C19 demonstrated low K(m) and high V(max) values for the metabolism of both model compounds, which supports their role as the primary enzymes that regulate metabolism at low-level human exposures to OPs. With K(m) and V(max) values of 0.61 microM, 4827 pmol/min/nmol P450 and 0.81 microM, 12,544 pmol/min/nmol for formation of paraoxon and chlorpyrifos-oxon, respectively, CYP2B6 favored the desulfation reaction. CYP2C19 activity favored dearylation with K(m) and V(max) values of 0.60 microM, 2338 pmol/min/nmol P450 and 1.63 microM, 13,128 pmol/min/nmol for formation of p-nitrophenol and 3,4,5-tricholorpyrindinol, respectively. P450-specific kinetic parameters for OP metabolism will be used with age-dependent hepatic P450 content to enhance PBPK/PD models so that OP exposures can be modeled to protect human health in different age groups.